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Abstract. It is shown that strong field approximation (SPA) can be formulated in 
a gauge invariant manner order by order of the expansion with no need for various 
partitions of the Hamiltonian in different gauges. 

PACS numbers: 32.80.Fb, 32.80.Rm 
Gauge-invariance of the transition probabilities describing interaction of atomic and 
molecular systems with intense electromagnetic fields is a long-standing problem. The 
point is that gauge-invariance can be apparently broken by strong field approximation 
(SFA), which is most commonly used for the semi-analytic calculations of the amplitudes 
[1]. For example, it has been shown recently [2] that ionization rates and energy 
distributions of strong-field photoionization differ substantially in the Keldysh-Faisal- 
Reiss theory [3l HI [5] when calculated in the length- and velocity gauges. Dependence 
of the SFA amplitude on the choice of gauge has raised a question which gauge is most 
appropriate for a description of a particular ionization process, depending on the atom, 
molecule or even on the initial state. For instance, it has been shown in [6] that in the 
case of ionization of negative ions with a ground state of odd parity the length-gauge 
SFA agrees better with a numerical solution of time dependent Schrodinger equation 
(TDSE). On the other hand, velocity-gauge calculations seem to reproduce quite well 
experimental results on photodetachment of negative fluorine ions whereas to get 
a similar agreement in the length-gauge would require alteration of measured peak 
intensity by a factor ca. 1.45. 

Whereas gauge invariance of the transition amplitudes (up to a phase factor) and 
probabilities cannot raise controversies as such, this fundamental property may be 
violated by various types of approximate calculations. It has been shown in [8] that the 
velocity- and length-gauge SFA amplitudes in the dipole approximation are equivalent in 
all orders of the expansion, provided appropriate initial- and final-state partitions of the 
Hamiltonian are chosen. These partitions in general differ from the "natural" partition 
into the unperturbed atomic Hamiltonian and the perturbation describing iteraction 
with the laser pulse, 

H{t) = H,t + F{t), (1) 
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where 

H.t = Ho + Kt (2) 

is the atomic Hamihonian. In Equation ([2]) Hq = p^/2 is the free Hamiltonian with p 
- the momentum operator, \4t is the binding (atomic) potential, and F(t) describes 
interaction with electromagnetic field in dipole approximation, with the property 
F(t) for t ±oo. For strong external fields it is rather the atomic potential 
which constitutes a perturbation, at least during action of a strong pulse, so that also 
another partition may be considered as "natural" 

H{t) = Hp{t) + K„ (3) 

where 

HF{t)=Ho + F{t). (4) 

In [8] Equation (jlj) has been chosen as the final-state partition, with the initial-state 
partition different from and given by (in atomic units) 

H{t) = H,{t)-A{t)-r, (5) 

where f is the coordinate operator, A{t) is the vector potential in dipole approximation 
and 

H,{t) = + + A{t) . r. (6) 

The momentum operator in Equation corresponds to the canonical momentum 
p = — iV. It has to be distinguished from the kinetic (gauge invariant) momentum 
Pkin = p + A(t). Both coincide at asymptotic times, since we assume that the 
vector potential is then zero. Wave function fulfilling Schrodinger equation with the 
Hamiltonian ([6]) has the form [9] 

$(r,t) =e-'"-^(*V(^,^), (7) 
where ip{r^t) is the atomic electron wave function obeying Schrodinger equation 
with ^at- 

A wider class of gauge transformations, of which the velocity- and length-gauge 
are particular cases, was considered in [in] 0. The main conclusion of [lOj is that 
constructing transition amplitudes describing interaction with strong external fields it 
is necessary to consider both the gauge and appropriate partition of the Hamiltonian. 
In general, expansions of the amplitude in different gauges can lead to the same results 
only in the limit of infinite series. 

In this letter we propose an alternative form of SEA expansion using a guiding 
principle that it is an expansion in powers of atomic potential, so that powers of V^t 

X The "traditional" minimal coupling velocity gauge obtained by replacement p — > p + is in [10] 
named as radiation gauge. In the "true" velocity-gauge the interaction Hamiltonian is p ■ A{t) without 
the A^{t) term. In this paper the velocity- gauge is understood "traditionally". 
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in every term in the series have to be carefully accounted for. The velocity-gauge 
Hamiltonian in the dipole approximation has the form 

Hy = ^[p + A{t)]^ + Kt = ^at + Fy{t), (8) 

where 

Fy{t)=p.A{t) + lA\t). (9) 

Transition to another gauge g is implemented by a unitary transformation of the form 
exp[ixg(r, t)]. Vector potential A, scalar potential (p and the state vectors \ip) transform 
according to 

Ag = A-Vxg, ^g = V + dtX9, 1^3) =e''^«|^), (10) 
and the transformed Hamiltonian reads 

dXgir,t) 



dt 



= \[p+ A{t) - VXgf + Kt - (11) 

Hamiltonian of interaction with the external field in the gauge g is 

-ix. _ 9xg{r,t) 



Fg{t) = e'^«i/ate-^^« - H^t + e'^^Fye 



dt 



1 1 f)^ 

= Fv(t) --ip-VXg + VXg ■ P) - VXg " A{t) + - (V^,)' - (12) 

For gauge transformations compatible with the dipole approximation Xg(^) ^) is at most 
linear in r. Whereas the vector potential should vanish at asymptotic times, it is not 
in general a necessary requirement for Xg) as can be seen, for instance, in the case of 
transformation from the "traditional" velocity-gauge (radiation gauge) to the "true" 
velocity gauge [TO] . 

Using partitions of the Hamiltonian defined by ([T]) and ([3]) one obtains two forms 
of integral equations for the time evolution operator JJ {t, t') 

Uit.f) = UAt,t')-i [' dtiUit,t,)F{t^)U^,{t,,t'), (13a) 

Jt' 

Uit,t') = UF{t,t')-i [' dtiU{t,ti)V,tUF{tut'), (136) 

Jt' 

where ?7at is generated by atomic Hamiltonian ([21) and Up - by the Hamiltonian (jl]). 
Transition amplitude from an initial bound electron state |0i(t')) at time t' to the final 
continuum state \(f>i{t)) (orthogonal to \(f)i)) at time t has the form 

M = (0f(t)|f/(t,t')|0i(O) = -i rrfti(0f(t)|f/(t,ti)F(ti)|0i(ti)), (14) 

Jt' 

where |0i(ti)) denotes the initial state propagated to the transient time ti by Ua.t(ti,t'), 
according to (113 ap . It has to be noted that due to partition (fl]), with H^t as the 
unperturbed Hamiltonian, the initial state is propagated to transient times always by 
the evolution operator ?7at, generated by = — + \4t- Gauge dependent objects 
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in expression (1141) are the complete evolution operator tJ{t,ti) and the interaction 
Hamiltonian F{ti). As can be seen from first part of Equation f|T^ transition amplitude 
is gauge independent (up to constant phase factor), since under the gauge transformation 

U{t, t') e'^«Wf/(t, t')e~'^^^''\ (15) 

For an asymptotic time t ^ oo the function Xg either vanishes or tends to a constant 
limit (independent of r), and vanishes for t' —oo [10]. For a short pulse with the 
vector potential vanishing beyond the time interval (ti, tf) it is sufficient to put t' and t 
as equal to t\ and tf respectively, since beyond this time interval the time evolution is 
generated by the unperturbed atomic Hamiltonian. 

As usual, transition amplitude is expanded in powers of the atomic potential by 
substituting consecutive iterations of fll36p into the second expression in (fl^ . First 
term in the SFA-expansion reads 

Mo = -i rrfti(0f^(t)|?7F(t,ti)F(ti)|0i(ti)), (16) 
Jt' 

where |(^f (t)) denotes SFA-approximated final state. Final state is usually approximated 
by a plane wave \q) P or by a plane wave distorted by asymptotic Coulomb phase [TT]. 
In the first case, which is more appropriate eg. for description of detachment from 
negative ions, U-p{ti,t)\(l)^{t)) is the Volkov state. Approximation of the final state by 
distorted plane wave should be used in the calculations of ionization of neutral atoms, 
when final electron moves in the long-range potential of residual ion. The initial and 
approximated final state are not orthogonal to each other. Using equation fulfilled by 

-i^UF{t,t,)=UF{t,t,)[Ho + F{t,)], (17) 
oti 

substituting the product f/pF from (fT7|) into (ITBi) and performing integration by parts 
one finds 

Mo = {^m\UfMm')) - (0^(t)|0i(t)) -i rrfti(0f-(t)|t/F(t,ti)Kt|0i(ti)), (18) 

Jt' 

where also the equation idt\4>i(t)) = [Ho + Va.t]\(pi(t)) fulfilled by the initial state was used 
at intermediate stage of the calculation. Calculation of the amplitude Mq in another 
gauge g, but still with the use of partitions ([1]) and ([3]) of the Hamiltonian, gives 

Mo^ = -2 rrfti(0^(t)|?7|(t,ti)F3(ti)|0i(ti)), (19) 
Jt' 

where Ul,{t,ti) = e'^^^*'>UFit,ti)e-'^^^*^'> . Substituting ^^(ti) from and using 

expression for Lfp one obtains from f|T9|) after simple though little tedious calculation 

-ie^>^«W f dt,{<pmUF{t,t,)e''^^^'^^VMti)) 
Jt' 

-ie'^^^W rrfti(0^(t)|f/F(t,ti)e-^^«(*^)Kt|0i(ti)). (20) 
Jt' 
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First two terms in (ITS!) and (I2U]) are equal provided that, apart of the property Xg(t') = 0, 
also Xgit) = for asymptotic times. This assumption narrows the class of gauge 
transformations considered in [10], but still leaves a wide class of transformations of 
the type Xg(r,t) = 7r ■ A{t), where 7 is an arbitrary real parameter. However, 
remaining contributions in (ITSl) and (|20l) . expressed as time integrals, are essentially 
gauge dependent, since Xgi't) 7^ for transient times. For strictly periodic long wave 
trains, such as considered eg. in [1], first two terms are zero distributions in the limit 
of asymptotic times, and what remains is a gauge dependent last term in ( |T8l) . Using 
the plane wave for 10^) and writing V"at as -f^at ~ -^0 one can then reproduce from (ITSll 
Reiss expression for the amplitude 

We shall show that consecutive terms in the SFA-expansion with the Hamiltonian 
partitions ([T]) and ([3]) are gauge independent if all contributions with given power of 
the atomic potential are groupped in one term in the expansion. Therefore, the term 
which is gauge dependent in (ITSl) . as linear in t4t5 will be absorbed by next term of the 
"conventional" SFA-expansion, given by 

Ml = {-if f'dt, f dt2{m)\MtM)%t%{t2M)F{ti)\Uti))- (21) 
Jt' Jtl 

Changing the order of integration, using again (fT7|) to express the product UpF and 
performing integration by parts we obtain 

Ml = i f dt2{<j)m\UF{t,t2)Ktmt2)) -i f' dt2{<pm\UF{t,t2)vMt2,t')mt')) 
Jt' Jt' 

+ (-i)2 f'dt2 r dt^{rmMtM)%tUY{t2M)%t\Hti))- (22) 
Jt' Jt' 

Adding Mq and Mi we see that the gauge dependent term in f|T8|) is cancelled by the 
first term of ( l22ll . and the last term in ( l22i) is already of the second order in the atomic 
potential and should be included into next term in the SFA-expansion. Thus 

Mo + Ml = M(°) + M(^) + 0(Kt), (23) 

where 

M(o) = {rMuAt,t')m')) - mt)\m). (24a) 

M« = -i f dtM{t)\lJ^{tM)%tMti,t')\Ut'))- (246) 

Jt' 

The terms M(°) and M^^) are separately gauge independent and their sum gives 
expansion of the transition amplitude up to and including contributions 0(Kt)- The 
first contribution M^^^ contains no atomic potential and describes direct ionization. The 
term M*^^) corresponds to the situation when emitted electron is scattered once by the 
atomic potential. It provides a simple picture of this process: initial state \(f)i{t')) at an 
initial time is propagated to the transient time ti under the influence of the pulse field, 
then "encounters" the atomic potential, and after that propagates to a final time again 
under the action of the pulse field only. 
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Next contribution to the amplitude can be obtained by substituting second iteration 
of fll36p into f|T^ . which gives 

M2 = {-if f dt^ f dt2 f dts{<j)t{t)\UF{t,ts)V,tUF{t3,t2)VMt^^ 

Jt' Jti Jt2 



Changing order of integration, using again flTTI) to express the product UpF and 
performing integration be parts we get for M2 

M2 = ^*rft3y[*'t^t2(0^(t)|t/F(t,t3)Ktf/F(t3,t2)Kt|0i(t2)) 

+ (-i)2 /*dt3 rdt2(0f(t)|f>F(t,t3)Kt?7F(t3,t2)Kt?7F(t2,Ol0i(O) 
Jt' Jt' 

+ (-i)3 fdts r dt2 rrftl(#(t)|?7F(t,t3)Ktf>F(t3,t2)Kt?7F(t2,tl)Kt|0i(tl)). (26) 
Jt' Jt' Jt' 

First term in fl26l) is equal to the OiV^^) contribution in fl22l) with opposite sign, and 
the last term in ( !26l) is of the order V^^. Therefore, up to and including contributions 

M = M^^^ + M^^^ + M^'^\ (27) 

where 

M(2) = (-i)2 f dt^ r dt2{4>m\MtM)V.tMhM)%tMt2MUt')). (28) 
Jt' Jt' 

and M*^°) and M^^-* are given by (124 al) and f l246[) respectively. It is easy to check that in 
an arbitrary order with respect to the atomic potential 

M(") = (-i)" ^* dtn fj dtn-i . . . fj dti 

X (0f^(t)|?7F(t,t„)Ktf/F(t„,tn-l)V;t . . . Ktf>F(tl,t')|0i(O)- (29) 

In the expansion of the transition amplitude 

oo 

M = M(°) + J2 M^^^ (30) 
fc=i 

every term is separately independent on the choice of gauge, where gauge invariance 
has been achieved by appropriate grouping of terms of the same order in the atomic 
potential. 

We have shown that the ^-matrix expansion including rescattering contributions 
to all orders can be formulated in a gauge invariant way provided that all terms of the 
same order in atomic potential are groupped together. In this approach every term of 
the expansion is separately gauge independent. A consequence of groupping together 
all terms with the same powers of atomic potential is that the first term in the SFA- 
expansion differs from the conventional one, since part of the latter had to be included 
into the 0(Kt) contribution. In every order of expansion gauge dependent parts from 
previous and next terms cancel. The new form of the expansion proposed here can be 
used to calculate transition amplitude in intense fields in the SFA-framework up to an 
arbitrary order in atomic potential with no necessity to choose a particular gauge. 
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